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ABSTRACT 

Observations have shown that continuous radio emission and also sporadic 
Ha and X-ray emission are prominent in singular, low-mass objects later than 
spectral class M. These activity signatures are interpreted as being caused by 
coupling of an ionised atmosphere to the stellar magnetic field. What remains a 
puzzle, however, is the mechanism by which such a cool atmosphere can produce 
the necessary level of ionisation. At these low temperatures, thermal gas pro- 
cesses are insufficient, but the formation of clouds sets in. Cloud particles can 
act as seeds for electron avalanches in streamers that ionise the ambient gas, and 
can lead to lightning and indirectly to magnetic field coupling, a combination of 
processes also expected for protoplanetary disks. However, the precondition is 
that the cloud particles are charged. 

We use results from Drift-Phoenix model atmospheres to investigate col- 
lisional processes that can lead to the ionisation of dust grains inside clouds. We 
show that ionisation by turbulence-induced dust-dust collisions is the most effi- 
cient kinetic process. The efficiency is highest in the inner cloud where particles 
grow quickly and, hence, the dust-to-gas ratio is high. Dust-dust collisions alone 
are not sufficient to improve the magnetic coupling of the atmosphere inside the 
cloud layers, but the charges supplied either on grains or within the gas phase as 
separated electrons can trigger secondary non-linear processes. Cosmic rays are 
likely to increase the global level of ionisation, but their influence decreases if a 
strong, large scale magnetic field is present as on Brown Dwarfs. 

We suggest that although thermal gas ionisation declines in objects across 
the fully-convective boundary, dust charging by collisional processes can play an 
important role in the lowest mass objects. The onset of atmospheric dust may 
therefore correlate with the anomalous X-ray and radio emission in atmospheres 
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that are cool, but charged more than expected by pure thermal ionisation. 
Subject headings: Brown Dwarfs, atmospheres, dust, ionisation, magnetic coupling 
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1. Introduction 

Both Brown Dwarfs and planets form clouds in their atmospheres. The radial extent of 
these clouds is mainly a function of gravity (Helling et al. 2010). High altitude clouds in the 
giant gas planet HD 189733b were inferred from transit measurements (Pont et al. 2008, 
Sing et al.2009). HD 189733b was the first exo-planet to be detected in polarized light which 
is likely to be caused by Rayleigh scattering on sub/im grains (Berdyugina et al. 2011). 
Lecavelier Des Etangs et al. (2008), too, suggest sub/zm dust grains as source of Rayleigh 
scattering in their analysis of the HD 189733b transit spectrum. The understanding that 
clouds are present in Brown Dwarf atmospheres was triggered by the disappearance of 
metal-oxide band (like TiO and VO) with decreasing effective temperature (see Kirkpatrick 
2005). The comparison with model atmospheres showed that cloudless models completely 
fail to explain the NIR colors and L-dwarf spectra (Chabrier et al. 2000, Burrows et al. 
2006, Saumond & Marley 2008). By today, various approximations of clouds as opacity 
source are in application (e.g. Burgasser at al. 2010, Currie et al. 2011) as the existence 
of clouds in Brown Dwarfs is well accepted. Photometric variability in L- and in T-dwarfs 
has been attributed to cloud evolution (e.g. Bailer- Jones 2008, Artigau et al. 2009), a 
finding which is supported by our work on turbulent dust formation (Helling et al. 2001, 
2004). Goldman et al. (2008), however, could not confirm spectroscopically the so-called 
'cloud-braking' across the L-T transition based on their model atmosphere interpolation. 

Scholz & Eisloffel (2004) argue that very low-mass objects (0.08 . . . O.25M ) show 
co-rotating, magnetically induced spots, which would indicate the presence of localised 
magnetic field concentrations. Objects at the boundary between M-dwarfs and Brown 
Dwarfs (0.08 . . . O.O9M ) were studied by several groups considering the activity change 
across the substellar border. The presence of X-ray and Ha emission (Berger et al. 2008, 
Reiners & Basri 2008) argues for a magnetically-heated corona, but the associated increase 
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in radio emission with decreasing mass suggests that the processes that are energising 
such a corona differ from those on higher-mass stellar objects (Hallinan et al. 2006, 2008; 
Antonova et al. 2008). In this case, the solar paradigm for coronal heating, which involves 
the build-up and release of magnetic stresses in the corona, may be ineffective. 

We have argued in Helling et al. (2010) that in the atmospheres of Brown Dwarfs and 
exoplanets, rapid ionisation events can occur, even in a cloud that is globally neutral and 
obeys dust-gas charge equilibrium. This can only happen, however, if the grains within the 
dust cloud are ionized. As one possible process, electron avalanche processes in the electric 
field of charged grains than lead to the formation of a streamer resulting in 10 13 . . . 10 14 
cm~ 3 free charges per initial electron (Dowds et al 2003, Li et al. 2007). Streamers are 
a growing ionisation front that propagates into non-ionised matter and are present in 
lightning ladders and sprites. 

The timescale on which these charges could recombine (the Coulomb recombination 
time scale) is longer than the time to establish a streamer over a large fraction of the cloud 
(Helling et al. 2010). Hence, a considerable degree of ionisation builds up in this time 
window during which the electrons could couple to the large-scale magnetic field. 

In this paper we evaluate the underlying assumption in this scenario - that the 
grains in substellar clouds are ionised. We concentrate on the influence of dust dynamics 
(sedimentation, turbulence) and kinematics (collisions) on the ionisation state of Brown 
Dwarf and giant planet atmospheres. Section [2] briefly summarises our kinetic dust cloud 
model and the results relevant for the topic of this paper. In Sees. [3j [3] we study the 
collisional energies of dust-gas and dust-dust collisions due to dust dynamics, and compare 
these results to the ionisation potential of solid compounds for an example of a brown 
dwarf and a gas planet atmosphere. Section [5] discusses potential lightning regimes in these 
atmospheres and Sect. [6] presents our conclusions. 
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2. Cloud formation model 

The formation of clouds is the formation of solid particles or liquid droplets in an 
atmosphere whose stratification is determined by gravity, radiation, and convection. Cloud 
formation in Brown Dwarfs and giant gas planets starts with the formation of seed particles 
(nucleation) . This differs from the cloud formation process on Earth (or solid planets in 
general) where dust grains that are carried into the atmosphere from the planet's surface 
act as condensation nuclei. In fully gaseous objects such as Brown Dwarfs or giant planets 
those nuclei need to be built within the atmosphere itself. Once these seeds have formed, 
many materials are thermally stable and grow a grain mantle by chemical surface reactions 
on the surface of each seed particle. These compounds form patches of different, adjacent 
materials on the grain surfaces as different solids can grow simultaneously. Hence, these 
clouds are made of grains that are composed of a mix of materials that depends on the 
local temperature and pressure (Fig. [TJ |2j top panels). As a consequence of these growth 
reactions, the mean particle size, (a) (Fig. [3j top panels), increases. 

Woitke & Helling (2003) and Helling & Woitke (2006) developed a model describing 
such heterogeneous dust formation by homogeneous nucleation, dirty growth, and 
evaporation. The formation of TiCVseeds onto which solid silicate and oxide compounds 
grow is considered. Woitke & Helling (2004) and Helling, Woitke & Thi (2008) applied 
this model to the formation of stationary clouds in oxygen-rich^ atmospheres including 
gravitational settling, element depletion, and convective element replenishment. The 
atmosphere model Drift-Phoenix (Dehn 2007, Helling et al. 2008a,b; Witte, Helling 
& Hauschildt 2009) couples this detailed kinetic model of dust cloud formation with a 

1 Oxygen-rich refers to a gas that contains more oxygen then carbon, hence the gas phase 
is dominated by a diversity of oxygen-binding molecules compared to a carbon-rich gas. The 
models used here apply the solar element abundances which are oxygen-rich. 
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radiative transfer code (Hauschildt & Baron 1999, Baron et al. 2003). We use the output 
of the Drift-Phoenix atmosphere simulations in this paper. 

In Drift-Phoenix the number of solids growing the mantle is restricted to seven 
(Ti0 2 [s], Al 2 3 [s], Fe[s], Si0 2 [s], MgO[s], MgSi0 3 [s], Mg 2 Si0 4 [s]) and to 32 surface 
reactions. Drift-Phoenix provides the local gas temperature T [K], the gas pressure p gas 
[dyn cm -2 ], the maximum convective velocity v™^ [cm s _1 ], and dust quantities such as 
the number density of dust particles n& [cm -3 ] of mean grain size (a) [cm] at each layer of 
the atmosphere. Because of the definition of the dust moments, pLj = V^ 3 f(V)dV (V 
- grain volume, V\ - volume of smallest possible grain), the dust moments Lj (j = 1, ... 4) 
that result from a Drift-Phoenix simulations are used to derive a representative grain 
size distribution function (for details see Appendix A in Helling, Woitke & Thi 2008). We 
consider in this paper a double-peaked grain size distribution f(a), 

f(a) = N 1 6(a-a 1 ) + N 2 6(a-a 2 ) (1) 



to calculate the relative dust velocities (Sect 4.1.2, 4.1.3) between grains of different 
sizes a in each layer of the atmospheric cloud (Fig. |3|. Parameters Ni, N 2 , «i and a 2 
are determined such that the resulting dust moments reproduce the solution of the dust 
moment equations from the Drift-Phoenix model atmosphere. Seed particles have the 
smallest size, a\, of all grains in the cloud and populate the uppermost cloud layers. As the 
cloud deck is defined as where the seed formation rate is maximum, small amount of very 
small particles form also above the cloud deck. 

While the dust grains are forming and increasing in size they fall under the influence 
of gravity into regions of increasing density (gravitational settling, sedimentation). The 
increasing gas density results in faster grain growth, such that the grains increase their size 
considerably over short distances, and in a decreasing sedimentation velocity. The relative 
velocities between dust and gas due to gravity can be considered as the equilibrium drift 
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velocity resulting from an equilibrium between gravity and friction with the surrounding 
gas (see Woitke & Helling 2003). The grains fall inward (i.e. downward) and they may 
reach a depth where the temperature is too high to sustain their thermal stability. For 
the example of a brown dwarf atmosphere, at T(p ga _ s 0.4bar) in Fig. [T] (top panel) the 
silicates evaporate, and instead, iron takes up w 75% of the grain volume. At even higher 
temperatures, i.e. higher gas pressures in Fig [I] (top panel), the dust evaporates completely 
which causes the decrease in grain size (Fig|3|. The same behaviour is found for the 
example of a giant planet atmosphere but at lower pressure of p gas ~ 0.0025bar . In both 
cases, the dust-to-gas ratios, pa/Pgas (black solid line, bottom Fig. [TJ [2]) is highest just 
above the evaporation zone where also most of the dust mass per cm 2 (red dotted line) is 
located. 

The dust influences the convective velocity v conv indirectly through the "backwarming" 
effect of its opacity on the local temperature gradient of the atmosphere. The backwarming 
changes the convective behaviour of the atmospheres because of a flatter or steeper 
local temperature gradient which then can suppress (5T/5r < 5T / 5r| a diabatic) or initiate 
(ST/5r > 5T/5r | a diabatic) convection locally resulting in detached convection layers. The 
observation of this effect is not new in model atmospheres and it can already be found 
in the first papers in Brown Dwarf atmosphere modelling (e.g. Tsuji 2002; Burrows, 
Sudarsky & Hubeny 2006). The convective energy transport in the Drift-Phoenix 
atmosphere simulations is considered in the mixing-length approximation from which we 
derive a large-scale convective velocity v conv . Helling (2005) showed that dust formation 
also influences the local velocity field by radiative cooling again due to the dust's large 
opacity. Hence, dust formation can sustain turbulence which enhances the relative velocities 
of dust grains on smaller scales than convection. However, the influence of the dust on the 
local turbulence spectrum is not taken into account in the Drift-Phoenix atmosphere 
simulation. 
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In this paper we use local quantities which have been self-consistently evaluated in the 
Drift-Phoenix atmosphere simulation (except for the turbulence spectrum) describing the 
formation of dust falling through a radiative and convective atmosphere. We investigate to 
what extent collisional processes are capable of ionising dust grains in clouds in atmospheres 
of Brown Dwarfs and giant gas planets. Figure [3] shows the dispersion in grain sizes a 
which produces a range of grain velocities at each height in the atmosphere (compare also 
Helling, Woitke & Thi 2006, Helling & Rietmeijer 2008). Various processes may influence 
the relative velocities between grains (Brownian motion, convection, gravitational settling, 
turbulence). We are particularly interested in whether collisional energies might then be 
large enough to overcome the ionisation potential of an electron on the grain surface, and 
hence if and how these processes can contribute to charge cloud particles in the atmospheres 
of Brown Dwarfs and planets. 

We note that the local (T, p gas )-structure can differ amongst model families due to 
differences in cloud modelling approaches (Helling et al. 2008). Even model atmosphere 
simulations for hotter objects differ enough to impact disk detection rates (Sinclair et al. 
2010). We are, however, confident that the Drift-Phoenix model atmospheres are a very 
reasonable attempt as the comparison to observations in the Dwarf Archieve demonstrates 
(Witte, Helling & Hauschildt 2011). 

We perform our study for a brown dwarf and a gas planet case. We choose two solar 
metallicity Drift-Phoenix models for a rather cool atmosphere of T c g = 1600K. The 
high-gravity case (logg=5.0) serves as an example for a brown dwarf atmosphere and the 
low-gravity case (logg=3.0) is the example for a gas planet atmosphere. This brown dwarf 
atmosphere model was also used in Helling et al. (2010). 
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3. Work functions for mixed-material dust 

The release of an electron from the gram's surface requires that it overcomes its 
ionisation potential, which is the energy needed to extract an electron from a surface and 
transport it to infinity. We assume this to be the work function in metals or the contact 
potential for insulators (see discussion in Poppe & Schrapler 2005). The values of the work 
function are mostly given for pure atomic materials in the literature (Desch & Cuzzi 2000, 
Kopnin et al. 2004) and range from « 2... 6eV (Fig. [5]) . However, our grains are made of 
well-mixed solid materials forming from an oxygen-rich gaseous environment, and they are 
only partially dielectric. Rosenberger (2001) states that the work function of multi-material 
grains is typically lower than the ionisation potential of pure atomic materials. We therefore 
conclude that the work function for grains of mixed materials is not well constrained, 
and we indicate the potential range of work functions in Fig. [5] (orange bar). The upper 
value may be somewhat overestimated leaving space for processes that we did take into 
account. Ideally, these electrons would be free to move around, but grains are observed to 
preferentially produce surface charges on their collisional partners rather than free electrons. 
We start our investigation with considering multi-electron emission due to dust collisions 
which may remain on the grain surface or escape if their energy is large enough. We see this 
process somewhat in analogy to secondary electron emission where more then one electron 
could be released during the collision with an electron or ion. During such collisions, 
backward scattering is more efficient in lattices of semi-conductor materials or insulators 
compared to metals which is of interest because our cloud model predicts particles made of 
a mix of materials. The secondary electron emission coefficient ranges from 2.4 for MgO[s] 
to 4.6 for Al 2 03[s] (NaClfs]: 6) and is highest for mixed materials (Ag-Cs20-Cs: 8; Niedrig 
1992, p359). We therefore adopt these secondary electron emission coefficients, i.e. the 
number of free electrons produced per collision, as guidance for mixed materials. We further 
note in analogy to the secondary electron emission that an increased collisional energy will 
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not necessarily cause a continuously increasing number of electron releases because the 
impact may affect deeper electrons in the solid which require a larger energy. 

4. Collisional ionisation 

We study ionisation events by collisions inside the dust cloud layers that form in Brown 
Dwarfs and planetary atmospheres. We investigate limiting cases in order to study if grains 
are charged in substellar clouds and hence, if they can act as seeds for other, more powerful 
mechanisms like e.g. electron avalanche processes. We compare the collisional energies to 
the ionisation energies of the dust grain surface. Three collision mechanisms are studied: 

- dust-gas collisions due to gravitational settling 

- dust-dust collisions by differential sedimentation of cloud particles 

- dust-dust collisions due to turbulence 

We discuss the effect of collisional processes only. During an elastic collision, the kinetic 
energy remains constant as the internal energy of the colliding bodies does not change. 
But the internal energy changes in the case of an inelastic collision e.g. due to chemical 
processes, and hence, the kinetic energy of the colliding grains must change. Such inelastic 
collision could cause fragmentation, erosion, sticking, deformation or evaporation (Guttler 
et al. 2010), besides the electron release discussed here. The occurrence of these processes, 
which are not considered in this paper, depends on the grain mass and relative velocity. 

4.1. Collisional energy and work function 

The collisional energy is calculated for two-body collisions which can either be 
dust-gas or dust-dust collisions. This collisional energy is then the kinetic energy for these 
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two-particle collisions, 



£col = T^redVrel- (2) 



Here, v rc i is the relative velocity between the collisional partners and m red is the reduced 
mass for a two particle collision. Thus for dust-gas collisions (dg): 

m d m g 

m Icddg — ■ , [6 

m d + m g 

and for dust-dust encounters (dd'): 

m d m d > 

m rcdA& = ■ (4) 

m d + m d > 

where m d = ^na 3 p d is the mass of a dust grain, m d i the mass of the colliding grain, and 
m g = p, the mean molecular weight in [amu] of the gas . 

Clearly, the collisional energy depends on the relative velocity of the collision. We 
therefore consider the above mentioned three possible types of collision separately in the 
following sections. 



4-1.1. Dust-gas collisions due to gravitational settling 

In Sec. 2 we have seen that, once the dust starts forming, each grain falls through 
its gaseous environment under the influence of gravity and friction. The equilibrium 
sedimentation velocity between gravitation and frictional force during sedimentation is 
reached quickly, and therefore, we can apply the equilibrium drift concept (Woitke & 
Helling 2003) to calculate the sedimentation velocity, v sed (a). The dust sedimentation 
velocity is relative to the gas phase and depends on the grain size a. We restrict this 
consideration to the large Knudsen number case (for a grain moving through gas) to be 
consistent with the quasi-static dust model applied in Witte, Helling & Hauschildt (2009) 
whose Drift-Phoenix model atmosphere results we will use. We, however, distinguish 
between the subsonic and the supersonic case where ct is the local sound speed: 
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< Ct 



(5) 




> ct ■ 



(6) 



The dust grains move with velocities higher than the sound speed in the very outer layers of 
the atmosphere due to the very low densities of the gas phase in these regions. A small-scale 
turbulent fluid field would influence the sedimentation velocity such that grains have more 
time and grow to larger sizes. 



Differential sedimentation due to different grain sizes provides a source of collisions 
amongst the dust particles themselves. Applying the above representation of the grain size 
distribution function (Eq. [T]) in each layer of the atmospheric cloud to Eq. [5] and Eq. [6] 
allows us to calculate relative sedimentation velocities, v^f (a), between grains of different 
sizes in each atmospheric layer caused by differential sedimentation. In order to illustrate 
the range of collisional velocities we calculate the relative velocities between the lower and 
upper grain size limits a% and a^- 



At high altitudes where p gas < 10~ 6 bar, only nucleation is possible. Therefore, the variation 
in grain sizes is zero and all grains have the same sedimentation velocity at a given height, 
hence v^^a) = 0. This behaviour is quite general as Witte, Helling & Hauschildt (2009) 
show in their Fig. 2 for a whole grid of substellar atmosphere models. The greatest 
relative velocity is achieved in the region of efficient growth where the grain size differences 



4-1.2. Dust-dust collisions due to differential sedimentation 




(7) 
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can be largest. However, the relative velocity drops with higher pressures because the 
Pgas-dependency dominates if a ^const because v sed (a) ~ a/p gas (or v sed (a) ~ \/a/p g as for 
v sed > c T ; Eqs.|||6f. 



4-1.3. Enhanced dust-dust collisions in a turbulent gas 

A further velocity component leading to inter-dust collisions results from a turbulent 
fluid field. The role of turbulence in the atmospheres of substellar objects such as Brown 
Dwarfs and planets has been studied in Helling et al. (2001). Convective mixing and 
gravity waves are mechanisms to drive a turbulent fluid field in a substellar atmosphere that 
otherwise would be damped by the viscosity of the atmospheric gas. Whereas the frictional 
coupling between the dust and gas is important in determining the terminal velocity during 
sedimentation, friction can also cause these particles to couple to the turbulent component 
of the gas velocity. Whether the dust grains couple to the turbulent motion of the gas 
depends on the relation between the frictional time scale, 77, and the turnover time, r t , of a 
turbulent eddy (e.g. Volk et al. 1980). Particles of a given size accelerate to their terminal 
velocity over a frictional timescale 77 (i.e., the grain's acceleration time scale r acc in Sect. 
2.5 in Woitke & Helling 2003). They can, therefore, couple to the velocity field of those 
turbulent eddies whose turnover time r* is larger than 77 (r t > 77). The frictional coupling 
time is for the subsonic and the supersonic case (Woitke & Helling 2003), 



3p gas c T 1 



2apd v scd 
3p gaa v 



jSt V— > Cr. 



For homogeneous, isotropic turbulence the energy dissipation rate e is constant for all scales 
(Kolmogoroff 's theory), and one finds from dimensional arguments that inside the inertial 
range 

e = Cxj, (9) 
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where u is the fluid velocity associated with a scale I. The constant C\ = 0.7 (Jimenez 
et al., 1993). We determine the locally constant energy dissipation rate e, for the 
largest scale of maximum kinetic energy assuming Z max = H p (r)/ 10 (pressure scale height 
H p (r) ~ 10 4 . . . 10 6 cm depending on radial distance r) and adopt the maximum vertical RMS 
velocity, u as proposed in Freytag et al. (2010, using their Eq. 5 with logV = logw max ), 

log W) = logv c m a nv + logr„ - (logp gas - log p^)/(H v /H p ) (10) 

with v™^v the maximum convective velocity and p™* the gas pressure at this maximum 



velocity taken from our Drift-Phoenix atmosphere simulation. Equation (10) 
parametrises the velocity of the large scale that drives the turbulent fluid field in the 
atmosphere according to the local thermodynamic conditions. Such large scale motions 
can be either convection or gravity waves according to Freytag et al. (2010). Freytag et 
al. (2010) derive parametrisations for H v /H p (velocity scale height of the wave amplitude; 
their Eq. 3) and r v (the ratio of maximum convection energy to wave amplitude = 'mixing 
efficiency', their Eq. 4) from their grid of 2D convection models for substellar atmospheres. 
In this parametrisation, H v /H p and logr„ are constant for a given set of stellar parameters. 
Hence, we find H v /H p = 2.855 and logr„ = —1.348 for our cases of T e g = 1600K (solar 



metallicity) from their Eqs. 3, 4. With v™*v = 4-5 • 10 3 cm/s and p™g X = 2.0 • 10 7 dyn/cm 
(=20 bar) for \og(g) = 5,0, and v™* = 1.7 • 10 4 cm/s and p™ a s x = 3.0 • 10 5 dyn/cm 2 for 



log(g) = 3,0, Eq. (10) reduces to 



logw max (r) = 4.862 S_Pgas ( Drown dwarf, log(g) = 5.0) (11) 



log w max (r) = 4.800 - (giant planet, \og(g) = 3.0) (12) 



2.855 
2.855 

Hence, the maximum vertical RMS velocity has no strong dependence on the surface 



gravity. The energy dissipation rate e (Eq. [9]) depends on the atmospheric pressure p gas , and 



it therefore changes globally but it is constant in each atmospheric layer Ah[ (see Sect. 4.3). 
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The eddy turn-over time, i.e. the correlation time of turbulent fluctuations, is 

/ /Z 2 \ 1/3 

7t = -~[-] , (13) 



u 

and it depends on the locally characteristic scale I (eddy size). Since r t determines the 
lifetime of an eddy, smaller eddies have a shorter lifetime than larger eddies for a constant 
energy dissipation rate e. Hence, if we know e and /, we can determine r t . 

Whether and how the turbulent gas motion influences the dust dynamics depends on 
the relation between Tf and r t . Under certain conditions the turbulent gas interacts with 
the dust grains such that the dust grains acquire a relative inter-grain velocity component 
among them. In order to quantify the influence of turbulence on the dust dynamics, we use 
the formalism deduced by Morfill (1985) to calculate the turbulence-induced relative (drift) 
velocity between the grains in the cloud, Av ind t : 

Av indit ~ 

W^[(i+^r i +(i+^)- 1 ) 

- 2( (i + W (1 + W '- (14 > 



Equation 14 simplifies for grains of similar sizes (t/i ~ 172) to 

1 + 



A Mt ^W)»^. (15) 



Tt 

-2\ ', 



Following Morfill (1985) we represent (5Vg)2, which is any systematic velocity component, 
by Eq. 10 or Eqs fU| / jl~2j To calculate r t , we consider a typical eddy scale of I = 1cm to 



represent small-scale end of the turbulence spectrum. This choice is guided by the critical 
length scale / cr ; for which Tf = r t (Fig. |i| top panels) i.e. when the momentum exchange 
between gas and dust is most efficient, and consequently the turbulence- induced inter-grain 
velocity component is highest. We can argue with Eq. [13] that the eddies have already 
decayed before the grain fictionally couple to the gas (t/ > r t ) for eddy sizes I < Z cri . The 
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grains equilibrate with the gas, hence they move together with the turbulent eddies only if 
Tf < Tt for which I > Z cri follows. Figure [4] demonstrates that the behaviour is rather diverse 
in the Brown Dwarf case for our choice of I = 1cm, pointing to a multi-scale problem 
that only very approximately can be described by one pre-chosen scale (see also Helling 
2005). In the planetary atmosphere case, our choice of / = 1cm suggests that all grains 
are frictionally coupled. If I is treated in more detail, the impact of turbulent motions on 
dust-dust collisions will be larger in the planetary atmosphere because £ cr j < 1cm. Note, 
however, that grains will decouple from the gas due to gravity on the largest scales in the 
atmosphere. 

Figure lil (bottom panels) shows the turbulence-induced relative (drift) velocity Av in d,t 



for collisions of grains of different sizes (dashed line, Eq. 14) and for grains of the same size 



(solid and dotted line, Eq. 15) over the whole atmospheric pressure range. The turbulence 
driving velocity, u mSQl (p gas ), is plotted in comparison. There is a maximum of Avi n d,t in both 
atmosphere cases, brown dwarf and gas planet, whereas Av ind t clearly decreases towards 
higher pressures. In these inner atmospheric regions, Tf decreases because the gas density 
increases, and therefore the dust grains can adjust more quickly to the gas motion. 



4.2. Energy release during dust-gas and dust-dust collisions 



The collisional energies according to Eqs. [2j |5j [6j and Eqs. [7], and Eq. 15 are shown in 



Fig. [5j together with the ionisation potential (work function) for a mix of solid materials 
as it is expected for a grain surface based on our dust formation model. We explore the 
results of our considerations for two solar metallicity Drift-Phoenix model atmospheres 
T e flf = 1600K for log(g) = 5.0 (brown dwarf) and log(g) = 3.0 (planet). 

Figure [5] shows that the collisional energies can change by 1-5 orders of magnitude 
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throughout the cloud layer of the brown dwarf and even more within the planetary 
atmosphere. This change is partly directly due to the variation in grain mass which in 
turn depends on the local thermodynamic conditions in the atmosphere, and partly due to 
changing local conditions, such as the relative grain velocity, temperature, and pressure. 

The collisional energies due to dust-gas collisions (not shown) are not large enough 
to provide the ionisation energy, and hence, this kind of collisions are negligible for 
grain ionisation. Dust-dust collisions between grains of different sizes due to differential 
sedimentation ( "rain out" , short dashed brown line) reach the ionisation potential for grains 
only at rather low gas pressures (p gas = 10~ 5 . . . lCT 4 bar) in the brown dwarf atmosphere. 
Comparing Fig. [5] with Fig. [3] shows that this interval is related to strong changes in grain 
sizes due to grain growth and grain evaporation. 

For turbulence-induced collisions between grains of similar sizes (green solid line, 
Fig. [5] ) produce the largest collisional energies. Most of the turbulence induced inter-grain 
collisions have collisional energies larger than the ionisation potential and are therefore 
most efficient in ionising grains inside the clouds in substellar atmospheres. 

Our example of a giant gas planet atmosphere demonstrates that dust-dust collisions 
are only sufficient to overcome the dust material's ionisation energies in the inner and 
densest part of the cloud, and only of turbulence amplifies the collisional process. This 
effect will amplify if the inter-grain velocities are larger due to a more efficient frictional 
coupling then assumed in this paper. However, the geometrical height of the cloud covered 
is larger than in a high-gravity brown dwarf atmosphere (Fig. [5]). 



-19- 



4.3. Electron pressure equivalent due to dust-dust collisions 

In the first paper of this series, we argued that charged grains can act like small 
capacitors and initiate electron avalanche processes that produce large numbers of free 
electrons. If many of these events superimpose and enough electrons leak out of the 
streamers, the degree of gas ionisation could be increased to such an extent that the cloud 
might couple to a large scale magnetic field. In this paper, we have demonstrated that 
the cloud particles can be charged in substellar and planetary atmospheres on the basis of 
turbulence induced inter-grain collisions alone. The next step is to quantify the number 
of electrons that potentially can be produced by these collisions in the whole cloud. This 
would also allows to investigate if inter-grain collisions have the potential to produce enough 
charges to allow a coupling with a large scale magnetic field in the case that these charges 
are released into the gas phase where they have a greater mobility. 

The energies produced by dust-dust collisions are larger than the ionisation energy 
by several orders of magnitude in most of the brown dwarf cloud volume or a fraction 
of it in gas planets if the cloud is turbulent. Such high energies which also exceed the 
thermal electron energy (black dotted line, Fig. [5]), suggest that electrons could escape 
from the grain surfaces into the gas. Whether the electrons affected by dust-dust collisions 
remain within the dust phase or whether they indeed escape from the grain surface 
depends on their kinetic energy gained during the collisions and on the electronegativity 
of the surrounding gas phase constituents. In the case of destructive grain processes, the 
potentially charge- carrying dust surface increases and charges could escape more easily 
into the gas as evaporation processes are faster. Many more mechanisms can influence 
this process (see introduction to Sect. [4j in addition to profound uncertainties in the 
microphysics of dust charging. 

Desch & Cuzzi (2000) consider frictional charge transfer between grains of different 
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size. These authors assume grains of homogeneous composition, i.e. insulator grains vs. 
metal grains, which is unlikely in oxygen-rich astrophysical gases, at least in the case 
of substellar atmospheres. Figures [TJ |] show that grains are composed of a mixture of 
metal and insulator materials over a large fraction of the cloud's extension. Hence, a 
charge-conservation model for homogeneous particles can not straight forwardly be applied 
to substellar atmospheres. Further, the collisional energies due to turbulence induced 
inter-grain velocities exceed the level of ionization in some regions by a few orders of 
magnitude. That implies that the number of electrons which are supposed to be exchanged 
between dust grains by triboelectric effects according to Desch & Cuzzi might eventually 
be released off the dust phase, having gained a certain amount of kinetic energy during 
the collision. We therefore approximate these possibly free electrons in the form of a 
non-thermal electron pressure to quantify the effect. 

We consider the frequency of collisions between dust particles of radius a and dust 
grains of another size a' (0). We denote the size distribution of these collisional partners as 
f(a') and write the collision frequency u co \ as, 

/■oo 

z/ co i = 7i a 2 Tcd v rcl (a, a') f(a') da' (16) 

J &\ 

with /(a + a') the reduced grain radius and a\ is the radius of the smallest grain 

possible (for details see Sect. 5 in Woitke & Helling 2003). Assuming a delta-function-like 



representation of the grain size distribution function f(a') = n' d (a')5(a — a'), Eq. 16 
simplifies to 

u co i = 7ra^ cd v re i(a, a')n' d (a). (17) 



2 This frequency will increase or decrease depending on the charge number of the grain 
itself in a distance smaller then the Debye length. Diver & Clark (1996) and Stark et al. 
(2006) showed that the grain's electric field will be stronger if the grain departs from spherical 
symmetry. 
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Here, n d [cm 3 ] is the number density of the collisional partner of the grain of size a. The 



relative velocity v re i(a, a') between grains of different sizes given by Eq. [7], Eqs. 14 and 15 
From Fig. [5] we conclude that only Eqs. 14 and 15 are relevant so far. 



We consider dust grains of size a passing through consecutive cloud fractions of 
thickness Ahi(T,p) (£; Ahi = geometrical cloud extension AH c \ oud , Fig. 1 in Helling et al. 
2010) colliding with other grains at a certain height in the atmosphere. Note that Ahi will 
increase with decreasing log(g) as the geometrical extension of the cloud, AH c \ oud , changes 
with surface gravity (Fig. 1 in Helling et al. 2010). 

The total number of potentially released electrons per unit volume in a large-scale 
velocity field within each Ahi is then given by 

f°° Ah- 

nZ = N-v coX f(a)—f-da. (18) 
J a,i v \ a ) 

N is the number of electrons produced per collision. In case of a process of comparable 
efficiency such as secondary electron emission, N ~ 8 elementary charges would be produced 
per collision with a grain made of mixed materials (Sect. [3]). For other processes, N can be 
larger by orders of magnitude. Also smaller Ah± might reveal other, small-scale velocities 



as leading terms compared to the large-scale velocity considered in Eqs. (18) potentially 
increasing the number of collisions n*°j. 



Using f(a) = n d (a)5(a), Eq. 18 simplifies to 



n% = N-v^n i (a) : £^. (19) 



v BLU (a 

_sed 



The sedimentation velocity v (a) of the impinging (primary) grains is approximated by 
the equilibrium drift velocity in Eqs. [5] and |6j This approximation seems well justified, 
however, since the time for particles to accelerate to their terminal fall speed v sed (a) is small 
(Woitke & Helling 2003). Even if grains are stopped by a collision, the grains will continue 
falling under the influence of gravity and they will achieve their terminal drift velocity 
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again quickly. This does not take turbulent mixing into account which would increase the 
sedimentation time. Large-scale convective replenishment is taken into account as part of 
the atmosphere simulation where it impacts the dust formation process and by that the 
gravitational settling. 

We assume that the dust and the gas (including thermal electrons) have the same 
temperature, T gas = Td us t- Woitke & Helling (2003) showed that the liberation of latent 
heat of condensation and the heating by friction is well balanced by radiative cooling and 
inelastic collisions with the gas particles for grain sizes to be expected in clouds of substellar 
atmospheres. We further assume that the electrons thermalise with the gas or dust, hence 
T e = T gas . However, T e = T gas may lead to an underestimation of the electron pressure p e 
as electrons tend to have higher kinetic energies than T gas in particular when originating 
from non-equilibrium processes like in streamers. For a first approximation, we retain this 
assumption and calculate an electron pressure equivalent, 

Pe,col = n co \kTg ?iS . (20) 



We have shown in Sect. ^2 that only turbulence-enhanced dust-dust collisions are 
capable of achieving appropriate relative grain-grain velocities to free electrons from the 
outer grain lattice (Fig. [5]). These grains have sizes between « 0.01/im and ~ 0.5/mi 
(Fig. [3]). They are made of a mix of silicates (Mg 2 Si04[s], MgSiOsfs], SiC^s]) which takes 
up ~80% of the grain volume with ^15% iron inclusions (Figs. [TJ |2j). 



Figure [6] compares the atmospheric, thermal electron pressure p e as result of the 
Drift-Phoenix atmosphere simulation (solid black line) with the electron pressure p e , C oi 
resulting from turbulence-enhanced dust-dust collisions producing N=8 charges (guided by 
the example of secondary electron production, dotted line), N=10 3 (dashed line, Desch & 
Cuzzi 2000), and N=10 6 . Hence, if the number of charges produced during the dust-dust 
collisions is comparable to the low efficiency of secondary electron production, the local 
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electron pressure will not be affected in a brown dwarf atmosphere. It, however, can 
resupply the surface with charges, a precondition need for electron avalanche process and 
streamers as suggested by Helling et al. (2010). Other mechanisms are needed to increase 
the number of electrons released from the grain lattice to allow p ejC oi 3> p e , g as- Helling et al 
(2010) demonstrated that a release of N=10 6 electrons per atmosphere layer would allow 
to couple 50% of the cloud volume to a large scale magnetic field. We discussed electron 
avalanche processes as attractive possibility to achieve such a high number of free charges. 
In Sect. [5j we demonstrate that in which atmospheric regimes superposition of avalanche 
processes during dust-dust collisions are likely. 

In gas planet atmospheres, the situation becomes more favourable for dust-dust 
collisions to increase the local degree of gas ionisation as the thermal gas pressure is lower 
than in brown dwarf atmospheres. Also the atmospheric volume affected is larger as the 
cloud has a larger geometrical extension than in brown dwarfs. 

5. Discussion 

5.1. Lightning vs. coronal discharge 

We have demonstrated so far that dust-dust collisions can be energetic enough to 
release electrons from the cloud particle's lattice. These electrons can be expelled into 
the gas phase or can remain on the grain surface. In the case of the electron remaining 
on the grain's surface, an electric field will build up which can be strong enough to cause 
an avalanche process that produces an exponentially increasing number of free electrons. 
As the recombination time back onto the grain surface is rather large, these electrons can 
exist for a time in the gas phase, and hence locally increase the degree of ionisation for a 
certain time (Figs. |6l [7]). But with which efficiency does this occur in the cloud and could 
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it lead to the occurrence of lightning? A superposition of avalanche-streamer processes will 
lead to more and more free electrons for a short time period which then may be defined 
as lightning. An estimate of the time scale, t = n d /?; sed , on which dust particles pass 
through such a previously formed electron cloud, and hence, potentially initiate another 
electron avalanche, shows that a certain fraction of the cloud is prone to lightning-like 
discharge events in brown dwarfs and in gas planets, namely where r str < t < T^^ m in 
Figs [7j The borders of this lightning region are not very clearly defined as turbulence will 
decrease the effective sedimentation velocity (blue solid vs. dashed lines). If dust particles 
pass an electron cloud slowly and the electrons can recombine, lightning is less likely and a 
less powerful coronal discharge-like behaviour on smaller scales should be expected as long 

aa 'str - 'recom- 

The lightning and the coronal discharge regimes will exist in both brown dwarfs and in 
gas planet atmospheres. Our results suggest a hierarchy of lightning and coronal discharges 
with the lightning occurring at lower pressures in the upper part of the cloud. However, the 
maximum of the collisional dust-dust energies (Fig. [HJ is more likely to be located in the 
coronal discharge regime in our planetary atmosphere example, while dust-dust collisions 
act across both regimes in the brown dwarf example. How would then the cloud particles be 
charged if dust-dust collisions are unfavourable inside a potential lightning regime? Cosmic 
ray ionisation maybe a possibility. 

Cosmic rays (CR) seem an important source of atmospheric ionisation in the solar 
system planets. Observation of Earth clouds, however, suggest that the actual charge 
production is not overly efficient but nevertheless important for coagulation processes. 
Nicoll & Harrison (2010) determine 17 . . . 150e per particle in cloud edges on Earth which 
can be directly related to ionisation by cosmic rays. The charging of the cloud particles is, 
however, not a direct result of the impact of the high energy CR particle on the cloud but 
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rather of the ion current that develops from the CR ionisation of the gas above the cloud. 
These ions attach to the cloud particles. The question is whether galactic cosmic rays could 
be a global source of ionisation for extrasolar low-mass objects. 

6. Conclusions 

Dust clouds are an integral part of the atmospheres of very low mass objects like 
Brown Dwarfs and planets. Clouds determine the local chemistry by element consumption 
and they influence the radiative and convective energy transport by their large opacity in 
the atmospheres. The aim of this paper is to demonstrate that dust grains in Brown Dwarf 
atmospheres can be charged, and to investigate whether the presence of dust in Brown 
Dwarf atmospheres can contribute to its overall ionization level, a necessary condition for 
magnetic coupling of the atmosphere. For this purpose, we focused on collisional processes 
of the dust phase to cause additional ionization of the atmosphere, an aspect that has not 
yet been considered in earlier research. We, however, acknowledge that a large variety of 
micro-physical processes can be involved into the ionisation of a mineral cloud which have 
not yet been taken into account in our model. We find that collisional energies can be 
high enough to ionize the dust phase over the whole extension of atmospheric clouds of a 
late type Brown Dwarf (T eff =1600, log(g) = 5) and over a large part in a giant planet's 
atmosphere (T e g=1600, log(g) = 3) if the influence of the turbulent gas motion on the dust 
grains is taken into account. We interpret our results such that collisional ionization of 
the dust grains on its own does not provide, in the first place, an ionization level which is 
sufficiently high for magnetic coupling. A collisionally charged dust phase could, however, 
trigger secondary non linear charging processes such as electron avalanches that lead to 
lightning: In the case that the charging process is rather a charge exchange between 
colliding dust grains, the larger grains will be charged positively and the smaller grains 
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would carry away the negative charges. This may lead to large scale charge separation 
due to differential sedimentation and also differential response to the turbulent gas drag of 
grains of different mass/size. 

We suggest that although thermal gas ionisation may decline in objects across the 
fully-convective boundary, dust ionisation may take over in the lowest mass objects. The 
onset of atmospheric dust formation, the cloud depth and its particle characteristics may 
therefore correlate with the anomalous X-ray and radio emission in atmospheres that are 
cool, but contain highly charged clouds. 

In the context of this paper, it may be surprising that only a fraction of Brown Dwarfs 
seem to show activity in form of X-ray or radio emission and not nearly as much is known for 
extrasolar planets. The observed intermittency of the X-rays might be interpreted as sign 
for an intermittent dust cloud distribution as suggested by us in earlier works (Helling et 
al. 2004). Depending on the objects parameter like mass, effective temperature, metallicity, 
age, the atmosphere does change and so do the conditions for electrification which we are 
only beginning to explore. Our first study does suggest that cloud electrification is more 
pronounced in giant planets compared to the more compact brown dwarf atmospheres. 
However, both brown dwarfs and giant planets are prone to lightning events inside their 
dust clouds. Systematic investigation of this are part of our future work. 
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Fig. 1. — Material composition in volume fractions V s /Vtot (V s - volume of a solid s, Vtot - total 
grain volume) of the dust cloud grains at a gas pressure p gas in the atmosphere (top). The bottom 
panel shows the dust-to-gas ratio Pd/Pgas (solid line) and the dust mass column density Md ust 
[gem -2 ] (dotted line). All quantities are results of the hydrostatic Drift-Phoenix atmosphere 
simulation for T g = 1600K, log(g)=5.0 and solar metallicity. The gas pressure p gas is therefore a 
measure for atmospheric altitude (compare e.g. top panels in Fig. 2 in Witte, Helling & Hauschildt 
2009). 
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Fig. 2. — Same like Fig. [I] but for an atmosphere simulation for T c ff = 1600K, log(g)=3.0 and 
solar metallicity. 
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Fig. 3. — Height-dependent grain size distributions as result from the Drift-Phoenix atmosphere 
simulation. Top panels: grain sizes a, Bottom panels: number of grains N. (oi, N±) and (02, iVjj) 
represent small and large particle ensemble in Eq. [TJ Both model atmospheres are for T e g = 1600K 
and solar metallicity. The top figure represents a Brown Dwarfs with log(g)=5.0, and the bottom 
figure a giant gas planet with log(g)=3.0. 




Fig. 4. — Top panel: Critical length scale for which Tf = t%. Bottom panel: Turbulence induced 



relative inter-grain velocity component v^t applying I = 1cm (Eq. 13). The local sound speed ct 
and the large-scale gas velocity, u max which drives the turbulence in the gas are over-plotted. Both 
model atmospheres are for T e g = 1600K and solar metallicity. The top figure represents a Brown 
Dwarfs with log(g)=5.0, and the bottom figure a giant gas planet with log(g)=3.0. 
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Fig. 5. — Collisional energies [eV] due to turbulent dust-dust encounters (green) and due to 
dust-dust encounters (brown) during gravitational settling. The electrons thermal kinetic energy, 
^ , kin( e_ )' i s overplotted as black dotted line. The range of plausible dust work function is shown 
as orange bar. Both model atmospheres are for T e g = 1600K and solar metallicity. The top 
figure represents a Brown Dwarfs with log(g)=5.0, and the bottom figure a giant gas planet 
with log(g)=3.0. 
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Fig. 6. — Electron pressure resulting from secondary electron production of N=8 (dotted), N=10 3 
(short dashed), and N=10 6 (long dashed) electrons per dust-dust collisions compared to the atmo- 
spheric, thermal electron pressure (solid line) for two Drift-Phoenix atmosphere (T e g = 1600K, 
log(g)=5.0 (top) / log(g)=3.0 (bottom), solar metallicity) . The horizontal arrow shows in which 
Pg as -range dust-dust collisions are energetic enough such that E co \ > E wor ^. 
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Fig. 7. — Time scale comparison of a streamer establishing (r a t r , black solid line), recombination 
onto the dust grain (t^J^ for q = lOe, red dotted line; Fig. 4 in Helling et al. 2010), and streamer 
superposition time scale t = n d 1 ^ 3 /i> sed (blue solid, blue dashed for 10 2 x v sed ) based on the mean 
grain sizes from the Drift-Phoenix atmospheres (T e g = 1600K, log(g)=5.0 / log(g)=3.0, solar 
metallicity). Multiple encounters of streamer electron clouds with the electric field of charged grains 
can happen in the 'lightning'-regime, but not in the 'coronal discharge'-regime. The extension of 
the 'lightning'-regime decreases if turbulence slows down the grain settling process (blue solid line). 



